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, Abstract 

We show that the Hausdorff distance for two sets of non-intersecting 

O' ' hne segments can be computed in parallel in 0(log'^ n) time using 0(n) 

processors in a CREW-PRAM computation model. We discuss how some 
parts of the sequential algorithm can be performed in parallel using pre- 
viously known parallel algorithms; and identify the so-far unsolved part 
^ of the problem for the parallel computation, which is the following: Given 

' ' two sets of x-monotone curve segments, red and blue, for each red segment 

^ find its extremal intersection points with the blue set, i.e. points with the 

^ minimal and maximal i-coordinate. Each segment set is assumed to be 

intersection free. For this intersection problem we describe a parallel al- 
\^ gorithm which completes the Hausdorff distance computation within the 

stated time and processor bounds. 

cn 

^ 1 Introduction 



Evaluating the similarity of two geometric shapes is an important problem 
in different fields of computer science including computer vision and pattern 
recognition. One of the most natural similarity measures is the Hausdorff 
distance which is defined for any two compact sets. The directed Hausdorff 
distance between two compact point sets P and Q is defined as dH{P,Q) = 
maxpgp miuggQ g), where d{p,q) denotes the Euclidean distance between 
the points p and q. The (undirected) Hausdorff distance Dh{P, Q) is defined as 
maximum of the two directed distances: Dh{P; Q) = max {d^f (P, Q),dH{Q, P)}- 
If we can compute the directed Hausdorff distance, we can clearly compute the 
undirected distance within the same asymptotic time bounds. Efficient sequen- 
tial algorithms are known for the Hausdorff distance computation for P and Q 
being discrete point sets, sets of non-intersecting straight line segments [5] and 
for algebraic parameterized curve segments [4j. 

*This research was performed in scope of the DFG-project "Parallel algorithms in compu- 
tational geometry with focus on shape matching" under the contract number AL 253/7-1. 
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In this paper we show that the HausdorfF distance for two sets of non- 
intersecting line segments can be computed efficiently in parallel within 0(log^ n) 
time using 0(n) processors in the CREW-PRAM computation model. 

Due to the current trend in hardware development, where the performance 
increase is achieved through additional computing units (processor kernels) in- 
stead of increase in CPU speed, parallel algorithms have gained new popularity 
in the algorithm development. Some graphic cards support general computa- 
tions on graphics hardware (GPGPU) which comprises up to several hundreds of 
parallel processing units, which even more motivates for development of parallel 
algorithms. 

In this paper we use the PRAM as model of computation model instead 
of choosing one of the currently available hardware platforms, since we want 
to examine principal parallelism of the problem instead of concentrating on 
specific technical details of some hardware. We plan an implementation on 
a GPGPU platform. Certainly, it is not feasible to implement the PRAM- 
algorithm presented here directly, but its major underlying ideas should be 
useful. 

We briefiy recapitulate here the key steps of the sequential algorithm for 
computing the directed Hausdorff distance between two sets P and Q of non- 
intersecting line segments from (1) Construct the Voronoi diagram VD{Q) 
of the set Q; (2) For each endpoint p of a segment in P find its closest segment in 
Q using VD{Q) and compute the distance from p to that segment; (3) Determine 
the so-called "critical points" on the edges of VD{Q)\ (4) For each critical point 
q compute the distance from q to its nearest segment in Q; (5) Return the 
maximal distance of endpoints and critical points. 

The authors show that the critical points, i.e., the points where the di- 
rected Hausdorff distance dniPjQ) can be attained, besides the endpoints of 
the segments in P, are the intersection points of P with the Voronoi edges of 
Q. Furthermore, they prove that for each edge of VD{Q) only the extreme 
intersection points, i.e., the first and the last intersection point along the curve 
segment, are critical points (s. Figure [T]), thus reducing the total number of 
critical points to 0{n). For a;- monotone curves the extreme intersection points 
are the points with the minimal and maximal x-coordinate. A non-a;-monotone 
parabola segment can be split into two a;-monotone segments at the point with 
the vertical tangent. Thus, in the following we can assume that all Voronoi 
edges are x-monotone. 

For all but one steps of the above algorithm there is either a straightforward 
parallel version or, due to some previous work, it is known how to compute them 
in parallel: A parallel algorithm for computing the Voronoi diagram of a set of 
line segments is given in [llj . It runs in 0(log^ n) time on a 0(n) processor 
CREW-PRAM and uses the divide-and-conquer technique. For a given planar 
subdivision with n vertices a point location data structure supporting O(logn)- 
time queries can be constructed in O(logn) time on an EREW-PRAM with 
0{n) processors [13]. Thus, step (2) can be performed efficiently in parallel. 
Parallelization of steps (4) and (5) is straightforward, since the distance for 
each of the 0{n) critical points to the nearest segment in Q can be computed 
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(a) Trademark images represented by (b) Critical points (c) Hausdorff dis- 
sets of line segments tance 

Figure 1: An example of two sets of line segments. The Hausdorff distance is 
attained at a critical point on a Voronoi edge. 

independently. 

The only open question in the parallel computation of the Hausdorff distance 
is to determine the critical points, i.e., intersection points between segments 
and Voronoi edges, efficiently. The sequential algorithm uses the plane sweep 
technique, in which the endpoints of the line and parabola segments are sorted 
by a;-coordinate and a vertical line is swept from left to right across the scene. 
Each time an intersection between an edge e of VD{Q) and a segment in P 
is detected the edge e is removed from scene, thus preventing computation of 
further intersection points of the same edge. Another sweep is preformed from 
right to left to detect the intersection points with the maximal x-coordinate. 

Clearly, the sweepline technique is inherently sequential, and we therefore 
need different tools to compute the critical points in parallel. For the related 
intersection detection problem: Given n line segments in the plane, determine 
if any two of them intersect - there exist efficient parallel algorithms. In [1] 
a CREW-PRAM algorithm is presented that runs in O(logri) time and uses 
O(nlogn) processors. This result is improved in ^5,, where intersection detection 
is performed in 0(log?T.) time on a CREW-PRAM with 0{n) processors. The 
latter algorithm has optimal processor-time product because of the 51(71 log n) 
sequential lower bound for this problem ilS; . 

For the intersection reporting problem, i.e., finding all pairwise intersections 
between line segments, a CREW-PRAM algorithm with running time 0(log^ n) 
on 0{n + k/ log n) processors is given in |10j , where k is the number of reported 
intersections. Detecting all intersections between two sets of line segments, 
each of which is intersection free, can be performed in O(logn) time using 
0{n + fc/logn) processors in a CREW-PRAM model [T2^. 

We call the intersection problem arising in the Hausdorff distance compu- 
tation, which also may be of independent interest, the first-last intersection 
problem (defined below). Since edges of the Voronoi diagram of a set of line 
segments can be line or parabola segments, we want our algorithm for this 
intersection problem to work for more general sets of curve segments: 

Definition. Two sets of curve segments A and B are called well-behaved if 

every segment in AU B is x-monotone; no two segments of the same set have 
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a common point except possibly common endpoints; any two segments from dif- 
ferent sets intersect at most twice; all intersections between any two segments 
can be computed in constant time; and for every segment we can compute in 
constant time for a given x-coordinate the corresponding y-coordinate. 

Observe that if we split the parabola segments of the Voronoi diagram at the 
points with vertical tangent, then the set P of line segments and the Voronoi 
edges of Q are two well-behaved sets. The problem of finding critical points for 
the Hausdorff distance can then be formulated as: 

Problem (First-Last Intersection Problem). Given two well-behaved sets A and 
B of curve segments in the plane, for each segment a G A find the intersection 
points of a with the segments from the set B with the smallest and with the 
largest x-coordinate. 

All above mentioned line segment intersection parallel algorithms utilize the 
segment tree data structure [5], which is also used in this paper and is described 
in section [2] There is no trivial modification of the mentioned algorithms for 
intersection detection or intersection reporting problems which yields an efficient 
algorithm for first-last intersection problem. Here we present an algorithm that 
solves this problem in 0(log^ n) time on a 0{n) processor CREW-PRAM and 
thus prove the following theorem: 

Theorem 1. Let A and B be two well-behaved sets of curve segments in the 
plane with \A\ -\- \B\ — n. The first-last intersection problem for A and B can 
be solved in O(logn) time on 0{n\ogn) processors using 0(n\ogn) storage in 
the CREW-PRAM model. Alternatively, the problem can be solved in 0(log^ n) 
time on 0{n) processors. 

Theorem [l] together with the above mentioned previous work completes the 
proof of Theorem [2j 

Theorem 2. Given two sets P and Q of n line segments, such that no two seg- 
ments of the same set intersect, except possibly at the endpoints, the Hausdorff 
distance Dh{P,Q) can be computed in O(log^n) time on 0{n) processors using 
0{n\ogn) storage in the GREW-PRAM model. 

2 Segment Tree and Interval Tree 

In this section we briefly describe two data structures used by our algorithm. 

Segment Tree. Let 5 be a set of n line segments. Sorting the 2n endpoints 
of the line segments by a;-coordinate and projecting them onto the x-axis results 
in at most 2n-\-l intervals. A segment tree for the set 5' is a complete balanced 
binary tree T with the following properties: Each of the 2n + l intervals is stored 
at a leaf of T in sorted order. An internal node t; of T stores the interval that 
is the union of the intervals of its children. 
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A segment s G S* covers a node v if the interval of v is completely contained 
in the projection of s onto the x-axis, but the parent interval of v is not. Every 
node V oiT stores a cover-list - the list of segments that cover v. Additionally, 
every node v stores a list of segments that have at least one endpoint in the 
interval oi v - the end-list. An example of a segment tree with cover-lists of the 
nodes is given in Figure |2(a)| 




(a) A set of four segments and the corre- (b) A set of seven intervals and the corresponding 
sponding segment tree T with cover-hsts of interval tree with unsortcd interval lists at the 
the nodes of T. nodes. 

Figure 2: Examples of a segment and an interval tree. 

Since every segment is contained in the cover-lists of at most two nodes of 
each level of T and in at most two end-lists per level, the size of T is 0(n log n). 
A segment tree with cover- and end-lists can be constructed in O(logn) time 
on a 0{n) processor CREW-PRAM (see e.g. [1]). 

Interval Tree. Let / be a set of n closed intervals on the real line. An 
interval tree Tj that stores / is defined recursively as follows: If / = 0, Tj is 
a leaf. Otherwise, the root node v of T/ stores a reference value and a list 
of the intervals of / that contain r„. The left (right) child of v is an interval 
tree for the intervals in /, whose right (left) endpoint is strictly less (greater) 
than r„. Typically, the reference value is chosen to be the median of the 
endpoints in /. This ensures that the height of Tj is O(logn). The intervals of 
a node of T/ are stored twice: in one list sorted by the first endpoint, and in the 
second, sorted by the second endpoint (see e.g. [9 for a detailed description). 
An example of an interval tree is given in Figure [2 (b)[ 

An interval tree for a set of n intervals uses 0{n) storage and can be build 
in 0(71 log ?i) time with a sequential algorithm. Using the interval tree we can 
report all intervals that contain a query point in 0(logn -f k) time, where k is 
the number of reported intervals. 

For the parallel construction we can sort the 2n endpoints of the intervals and 
build a complete balanced binary search tree Tj on the values of the endpoints. 
Now, using the values in the nodes of the search tree as the reference values, 
each interval of / can independently find the highest node in T/ whose reference 
value it contains and assign itself to that node. Finally, we sort the node entries 
lexicographically twice: by (u, ii) and by (u, 12), where v is an identifier of a tree 
node and ii, i2 are the endpoints of an interval assigned to v. The intervals can 
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then be written to the two hsts of their corresponding nodes: in one sorted by 
start- and the other sorted by endpoint. This gives us a tree with the properties 
of the interval tree and O(logn) height. Thus, using a fast sorting algorithm we 
have: 

Lemma 3. For a given set ofn intervals an interval tree can he constructed in 
O(logn) time on 0{n) processors in the CREW-PRAM model. 

3 A Parallel Algorithm for the First-Last-Inter- 
section Problem 

Let A and B be two well-behaved sets (red and blue resp.) of curve segments in 
the plane with \ A\ + \B\ = n. Here we describe how to find for each segment a € 
A the intersection point with B with the minimal x-coordinate, the intersections 
with the maximal x-coordinate can be determined symmetrically. 

Our algorithm begins with the construction of a segment tree T for the set 
A\JB: 



Step 1. Build a segment tree T for AiJ B. For each node v of T construct 
separate cover-lists Ca{'v),Cb{v) for sets A and B respectively, and separate 
end-lists Ea{v), Eb{v). Sort Ca{v) and Cb{v) by y-coordinate for all nodes v 
in parallel. 



Since the initial sets are intersection free, the y-order within the cover-lists 
of each color is well-defined. 

Chazelle showed in that if two line segments of a set intersect, then in 
the corresponding segment tree there must be a node v such that either both 
segments are in C{v) or one is in C{v) and the other in E{v). In our example 
in Figure 2(a) the intersection of the segments S2 and S3 is of the first type and 
the intersection between si and S2 of the second. It is easy to see that a similar 
statement holds for well-behaved curves as defined in Section [T] Consider an 
intersection point p of curve segments s and s' and the leaf v of the segment 
tree that contains the x-coordinate of p. Then each of the two segments s and 
s' appears in the cover-list of one node on the path from v to the root of the 
segment tree. Let v' denote the highest, i.e., closest to the root, such node, and 
let w.l.o.g. s' be in the cover-list of v' . Then, either s is also in the cover-list of 
v' , or, since s covers a sub-interval (the interval of v) of v' but not the complete 
interval of v' or any of its ancestors, s must be in the end-list of v' . The same 
argument holds for every intersection point of s and s' and the corresponding 
leaf of the segment tree. 

In the two-set setting this means that if a red segment a and a blue segment 
b intersect, then there must be a node v in T such that either (1) a G Ca{v) and 
b e Cb{v), (2) a € Ea{v) and b € Cb{v), or (3) a G CAiv) and b e EBiv). The 
following steps of the algorithm deal with each of these three cases. Whereas the 
handling of the first two cases (Step 2) is a modification of the corresponding 
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steps of the algorithm in T?, the third case demands additional processing 
(Step 3). 



Step 2. For each node v of T and each segment a G Ca(v) U Ea{v) do in 
parallel: Find the neighbors of a in Cb{v) with respect to the y -order at the 
x-coordinate of the leftmost point of a within the interval of v. Compute the 
intersections of a with its neighbors, if any exist, and record the one with the 
minimum x-coordinate. 



Since we want to find the intersection point of a with the minimum x- 
coordinate and all curve segments are a;-monotone, we do not need to find all 
intersections of a within the interval of v but only those with the neighbors at 
the left border of the interval, or at the x-coordinate of the left endpoint of a, 
respectively. 

For type (3) of intersections we observe that for a given segment b G Eb{v) 
we can easily determine the lowest, ai, and the highest, a2, red segment in Ca{v) 
intersected by 6, such that the intersection points have x-coordinates inside the 
interval of w, using binary search on Ca{v)- Clearly, all red segments in Ca{v) 
between oi and 02 are also intersected by h forming a set of consecutive ranks 
in Ca{v) with respect to its ascending order in y-direction. This set we call the 
rank interval of 6, which has a constant size representation by 01,02. 

In the following we are going to find for each node u of T the set I{v) of the 
rank intervals for all b £ Eb{v). Then we process and narrow the rank intervals 
in I(v) with the purpose to include at most O(logn) intersection points for each 
a e Ca{v). 

The purpose of the next step is to avoid multiple intersections of a blue 
segment with a red one in its rank interval. Observe that there are three pos- 
sibilities for the number of intersection points of b with its lowest and highest 
intersected red segments oi and 02: (1) Both red segments are intersected once 



(Fig. 3(a)). In this case every red segment within the rank interval is also in- 
tersected exactly once. (2) One red segment is intersected once, and the other 
- twice (Fig. |3(b)[ ). Let oi have one intersection point and 02 two. Then all 
red segments intersected by b once are the consecutive red segments following 
ai, all red segments with two intersections are consecutive segments preceding 
02, and all twice intersected segments yield nested intervals on b. Thus, when 
we split b as in Step 3.1(c), one of the new segments has the same rank interval 
as b and for the other we need to determine the new lower end of the interval 
(the upper end is 02). (3) Both red segments are intersected twice. Then the 



x-intervals of their intersection points are either nested (Fig. |3(c~ ) or disjoint 



(Fig. |3(d)[ ) . If the intervals are nested, then all red segments of the rank interval 
are intersected twice in same manner, and both new subsegments of b have the 
same rank interval. Otherwise, all red segments of the rank interval have to 
intersect the middle subsegment of b but not necessarily the end subsegments. 

Step 3.1. For each node v in T and each b € Eb(v) do in parallel: 
1. Find the rank interval of Ca{v) intersected by b. 
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(a) One intersection (b) One intersection on (c) Two intersections (d) Two intersections 
on both ends one end and two on the on both ends - nested on both ends - dis- 

other joint 



Figure 3: Possible configurations of the intersection points of b with the first 
and last red segments of its rank interval. 



2. Compute the intersection points of b with the lowest and with the highest 
intersected segment in Ca{v). 

3. If one interval end yields one and the other two intersection points (Fig. \3(b)\ l, 
split b between the intersection points of the latter rank and find the new 
rank interval for one of the resulting segments (see below). 

4- If both interval ends yield two intersection points and the intersection in- 
tervals are: 



(a) nested (Fig. 3(c)^ - split b between the innermost intersection points, 



both resulting segments have the same rank intervals as b. 

(b) disjoint (Fig. \3( - split b between each pair of intersection points. 
The middle segment has the rank interval of b the intervals of the two 
end segments have to be determined. 

5. Record the resulting interval(s) in the interval set I{v). 



Now we can be sure that each blue segment intersects each red segment in 
Cyi(f) at most once within the slab of v. 

Step 3.2. For each node v G T in parallel construct an interval tree Tj{v) for 
the set of rank intervals I{v). 

The blue segments whose rank intervals are stored in the same node u of 
Ti{v) all intersect the same red segment a e Ca{v) - the segment with the 
reference rank of m, hereafter called the reference segment of u. Thus, we can 
order them by the a:-coordinate of their intersection with a. If two blue segments 
bi , 62 intersect some red segments the order of the intersection points will be 
the same on all of these red segments. Therefore, if 61 's intersection with a 
has a lower x-coordinate than that of 62 , we can remove from the rank interval 
of 62 those elements that are in the interval of bi without losing significant 
intersection points. 

Step 3.3. For each node v in T and each node u in Ti{v) in parallel do: 
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1. Sort the rank intervals of u by the x-coordinate of the intersection with 
the reference segment of u. 

2. For the resulting ordered sequence compute the prefix-maxima of the high- 
est ranks and the prefix-minima of the lowest ranks. 

3. Replace the i-th interval [r^^i, ri_2\, i = 2, 3, . . . , by two intervals 

[r,;^i, niin,;_i — 1] and [maxi_i +1, ?'i,2]; where mhvi^i = minfc=i_...j_i rfc^i 
and niaxi_i = max^^i i_i 2 (see Fig. 

4-. Discard all empty intervals, i.e., intervals [ri,r2] where r2 < ri. 

The prefix-minimum (min^) and the prefix-maximum (maxj) for the i-th 
interval in the ordered sequence from Step 3.3 give us the range of ranks of the 
red segments that are already intersected by at least one of the corresponding 
blue segments hi, ... ,bi before they (possibly) encounter an intersection with the 
segment 6^+1. Thus, we can remove the [mini,maxi] interval from the interval 
of bi+i. 

After the interval spliting and shrinking Tj(v) is no longer an interval tree, 
but the remaining intervals have the following properties: The rank intervals 
stored in one node of Ti{v) are now disjoint, i.e., the rank of every segment 
a G Ca{v) is contained in at most one interval of a single node. The number 
of intervals in Tj{v) is at most twice the original number. Since every rank is 
stored in at most one node of each level of Tj{v), the rank of every a S Ca{v) 
is contained in at most O(logn) intervals in Ti{v). An example of such interval 
reduction for one node of the interval tree Tj (v) is given in Figure El 




fes [2,7] 1 7 [7,7] 



Figure 4: Segments in Ca{v) and Eb{v) corresponding to one node of Ti{v) with 
the reference rank 4. In the right picture the removed parts of the intervals are 
denoted by dashed lines. 

Finally, we can reorganize Ti{v) and compute the intersection points: 

Step 3.4. For each node v ^ T in parallel rebuild the interval tree Ti{v) for 
the new intervals. 

Step 3.5. For each red segment a £ Ca{v) in parallel find all b G Eb{v) whose 
intervals in Ti{v) contain the rank of a. Record the intersection point with the 
minimal x-coordinate. 
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The last step is to select the intersection points to report from the candidate 
intersections: 

Step 4. Sort the candidate intersection points computed in Steps 2 and 3 lexi- 
cographically by (oi, X- coordinate) , where Ui is an identifier of the red segment in 
A. For each a G A report the intersection point with the minimal x-coordinate. 

In the next section we prove that the described algorithm correctly reports 
the first intersections of each segment in A and do the counting to show the 
workload claimed in Theorem [TJ 

4 Proof of Theorem [I] 

The correctness of the algorithm is based on the fact that in Steps 2 and 3 
of the algorithm we consider all types of possible intersections and for each 
segment a G A discard only such intersections of each type which cannot have 
the minimal x-coordinate. We elaborate on the latter claim: 

Consider the intersections of type (1) and (2), i.e., intersections within the 
vertical slab of a node w of T of a segment a in CAiv) or in Ea{v) with the 
segments in Cb(?;), which are handled in Step 2 of the algorithm. Since there 
are no blue-blue intersections, the segments of Cb{v) partition the vertical slab 
corresponding to the x-interval of?; into (curvilinear) quadrilaterals. By locating 
the neighbors of a in Cb{v) at the leftmost point of a within the vertical slab 
of V we find the quadrilateral in which a "enters" the slab of v. Clearly, the 
first intersection for a within the slab of v (if any exists) is with one of its 
blue neighbors. Because all curve segments are a;-monotone, this is also the 
intersection with the lowest x-coordinate within the interval of v. We can safely 
discard all further intersections of a with Cb [v) ■ 

In Step 3 we initially capture all intersections of type (3) in the rank intervals 
without explicitly computing each intersection. Since our curve segments are 
well-behaved and we split the blue segments into subsegments that intersect 
each red curve at most once, if two blue sub-segments intersect both two red 
segments oi and a2, then the x-order of their intersection points is the same for 
ai and 02. Thus, when we exclude the rank of a red segment a from the interval 
of a blue sub-segment b in Step 3.3, we know that there is another blue sub- 
segment, whose intersection with a has a lower x-coordinate. The intersection 
of a and b can then safely be discarded. 

Next we prove the running time and the resource usage of the algorithm 
claimed in Theorem [1] Step 1 can be performed in O(logn) time on 0{n) 
processor CREW-PRAM, see [1]. A segment tree uses O(nlogn) space. 

In Step 2 we assign to each segment a G A a, processor which traverses 
the segment tree T and for each node v, with a e C^(f) U Ea(v) performs 
a binary search on Cb{v) to determine the neighbors of a. The intersection 
computation is performed in constant time. This gives us O(logn) time for 
each of the O(logri) levels of T, i.e., O(log^n) time with 0{n) processors for 
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all candidate intersections of type (1) and (2). Alternatively, we could assign a 
processor to each occurrence of a segment a in cover- and end-lists of T and find 
all candidates of type (1) and (2) in O(logn) time using 0{n\ogn) processors. 
The first variant gives us 0{n) candidate intersections in Step 2, the second - 
O(nlogn). 

In Step 3.1 using one processor for each b £ B we can find in time 0(log^ n) 
all rank intervals for the intersections of type (3) in T. Since each b occurs in the 
end lists of at most two nodes per level of the segment tree T, each b produces 
O(logn) such intervals. Sorting all 0(n log n) boundaries of the rank intervals 
lexicographically by (w, r), where v is the node of T that produced the interval, 
and r is a boundary of the interval, we get the grouping of the intervals by the 
nodes of T and sorting the interval boundaries in one step. This sorting can 
be performed in O(log^n) time using 0{n) processors. Within the same time 
and processor bounds we can construct all interval trees in parallel (Step 3.2). 
Since an interval tree uses linear storage, the total space requirements remain 
O(nlogn). 

The operations in Step 3.3 involve sorting lists with a total number of 
O(nlogn) elements, and performing prefix-max and prefix-min computations 
on these lists. These operations stay within the stated resource bounds. Split- 
ting of the intervals and eventually deleting some of them can be performed 
independently for each interval, i.e., in O(logn) time by 0{n) processors in to- 
tal. The number of the intervals is at most doubled in Step 3.3, thus the storage 
requirement remains 0(?ilogn) and Step 3.4 has the same resource requirements 
as Step 3.2. 

In Step 3.5 using one processor for each red segment a we can traverse the 
segment tree T and in each node v, such that a € Ca{v) find all intervals in Tj{v) 
that contain the rank of a in Ca{v). As mentioned above, there are at most 
O(logn) such intervals in Ti{v). Each interval corresponds to a blue segment in 
Eb{v). Instead of reporting all intervals we compute the intersection of a with 
the corresponding blue segment and keep track of the intersection point with 
minimal x-coordinate. The computation of all intersection points of a stored 
in T/(w) is performed in O(logrt) time, which gives us total of O(log^n) time 
over all levels of T. The number of candidate intersections produced in Step 3 
is 0{n) - at most one for each red segment. 

Alternatively, all operations of Step 3 can be performed in O(logn) time on 
O(nlogn) processors. 

In Step 4 we sort 0{n) candidate intersection points by (a,a;), where a is 
the number of the red segment and x is the x-coordinate of the intersection 
point, and take the minimum for each red segment. Step 4 can be performed in 
O(logn) time using 0{n) processors. Which concludes the proof of Theorem[l] 

Throughout the algorithm we clearly use concurrent read operations, but 
for each write operation for each processor we can independently determine 
a unique storage slot. Thus the algorithm is designed for the CREW-PRAM 
model. 
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5 Concluding Remarks 



The algorithm presented here can further be accelerated to perform Step 2 
in O(logn) time using 0{n) processors by applying the fractional cascading 
technique [8], [5]. This technique simplifies the iterative search for a key in 
multiple ordered lists and allows to perform the search of a segment a in all 
Cs-lists of T in 0{\ogn) time, i.e., taking 0{\ogn) time on 0{n) processors for 
Step 2. Steps 1 and 4 have already these time and processor bounds. But it is 
not clear whether the performance of Step 3 can be improved. 

A further interesting related problem is matching geometric shapes under 
transformations (e.g., translations, rotations, scalings) with respect to the Haus- 
dorff distance, i.e., find a transformation of one of the shapes such that the 
Hausdorff distance is minimized, for sequential algorithms see [2]. 

Another interesting distance measure is the Frechet distance [3]. A paral- 
lel solution to the decision problem, i.e., testing whether the Frechet distance 
between two polygonal curves is at most some given value e, should be uncom- 
plicated, since the sequential algorithm uses divide-and-conquer technique. We 
are currently working on a parallel algorithm for the computation problem, for 
which the sequential algorithm uses parametric search. 
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